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(~| I Abstract 



In the Next-to-Minimal Supersymmetric Standard Model (NMSSM), the presence 
^ I of light pseudoscalars can have a dramatic effect on the decays of the Standard Model- 

like Higgs boson. These pseudoscalars are naturally light if supersymmetry breaking 
f\| . preserves an approximate U{1)r symmetry, spontaneously broken when the Higgs 

Cn I bosons take on their expectation values. We investigate two classes of theories that 

t^ ' possess such an approximate U{1)r at the mediation scale: deformations of gauge 

and gaugino mediation. In the models we consider, we find two disjoint classes of 
f^ i phenomenologically allowed parameter regions. One of these regions corresponds to a 

limit where the singlet of the NMSSM largely decouples. The other can give rise to a 
Standard Model-like Higgs boson with a dominant branching into light pseudoscalars. 
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1 Introduction 

The Minimal Supersymmetric Standard Model (MSSM) is a leading theory for new physics 
at the weak scale, providing for a stabilization of the weak hierarchy, the unification of 
gauge couplings, and an excellent dark matter candidate [1]. In contrast to the Standard 
Model (SM), the MSSM possesses a pair of electroweak Higgs doublets. These give rise to five 
physical degrees of freedom instead of the single SM Higgs boson. Despite this complication, 
throughout much of the allowed MSSM parameter space one of the MSSM Higgs bosons 
behaves in very much the same way as the SM Higgs [2]. 

The MSSM is not without its puzzles. In particular, it possesses a dimensionful parameter 
in the superpotential that marries the two Higgs multiplets together: W D ^HuHd. This 
parameter must be of order the weak scale in order to achieve proper electroweak symmetry 
breaking. One approach to this problem is the Next-to-Minimal Supersymmetric Standard 
Model (NMSSM), which consists of the MSSM augmented by a singlet chiral superfield [3lll]. 
This singlet superfield leads to two new physical Higgs states: a neutral scalar and a neutral 
pseudoscalar. The appearance of these states, along with their mixing with the MSSM 
Higgs states, can significantly modify the Higgs phenomenology at the Large Hadron Collider 
(LHC) and other coUiders [11 E El E E] . 

The LHC Higgs boson signatures of the NMSSM are particularly different from the SM 
(and the MSSM) when the light SM-like Higgs boson can decay into pairs of very light 
mostly singlet Higgs pseudoscalars. (For a recent review of these decays see Ref. [9j.) It 
was observed in Ref. [10] that if this decay mode dominates over the standard SM modes 
such as /i ^ 66, the LEP-II bound on the SM-like Higgs mass is lowered below 114 GeV. 
Masses of 110 GeV or below are possible if the pseudoscalars decay primarily into bottom 
quarks [9l |1T], and Higgs bosons as light as 90 GeV can be consistent with the LEP data 
if the pseudoscalars decay primarily into tau leptons. Since the dominant contribution to 
fine-tuning within the MSSM comes from the Higgs boson mass bound, reducing the bound 
on the SM-like Higgs mass in this way can conceivably ameliorate the MSSM Higgs-sector 
fine-tuning problem. Recent discussions of this point have appeared in Refs. [TOl 



Often, the Higgs pseudoscalars of the NMSSM are too heavy for the SM-like Higgs to 
decay into them. An important exception occurs when the theory has an approximate 
continuous symmetry under which at least one of the Higgs states is charged. In this case, a 
light pseudoscalar arises as the pseudo-Nambu-Goldstone boson (pNGB) of the approximate 
symmetry when it is broken in the course of electroweak symmetry breaking. A promising 
candidate for this symmetry is a U{1)r [T3l [HI [15], which is exact in the limit where the 
trilinear A terms of the Higgs sector vanish. 

This symmetry is realized approximately in a natural way if the mediation of supersym- 
metry breaking is dominated by gauge interactions, and the Higgs sector A terms vanish at 
the mediation scale. In this case the dominant source of U{1)r violation comes from the 
(Majorana) gaugino masses. This breaking is communicated to the Higgs fields in the course 
of renormalization group (RG) running between the mediation scale and the electroweak 
scale. 



Two ways to mediate supersyminetry breaking through gauge interactions with vanishing 
A terms are gauge mediation (GMSB) [16], [TTl [H] and gaugino mediation (gMSB) p^ [20]. 
In both of these mechanisms of SUSY breaking, it is challenging to generate both the fi 
and the B/i terms with the correct relative size in the MSSM. One way to address this 
challenge is to add a singlet superfield whose vacuum expectation value (VEV) generates 
these terms, as is done in the NMSSM. Unfortunately, previous work on combining gauge or 
gaugino mediation with the NMSSM has found that the effective shielding of the singlet sector 
from the supersymmetry breaking also makes it difficult to obtain an acceptable pattern of 
electroweak symmetry breaking [21]. A number of works have investigated extensions of 
minimal gauge mediation in the hope of relieving this tension [22l [23| [2] 



In the present work, we study the LHC Higgs boson phenomenology of two simple 
extensions of minimal gauge and gaugino mediation within the NMSSM. The first extension 
consists of relaxing some of the GMSB boundary conditions on the singlet-sector soft terms. 
This might arise from a direct coupling of the singlet sector to the messenger states, or if 
there is additional source of supersymmetry breaking coupled solely to the singlet [26] . Such 
a modification of gauge mediation was considered recently in Ref. |27j while this work was in 
preparation. We extend and expand upon their results within gauge mediation, and apply 
the same modification to gaugino mediation. The second extension we consider involves 
adding new charged vector-like states to the theory [22l [21] . These modify the low-energy 
soft parameters of the NMSSM through their effects on the renormalization group (RG) 
running. Before discussing these modifications in detail, we introduce notation and review 
why a modification of the minimal gauge mediation scenario is needed. 



1.1 The NMSSM Higgs Sector 

The Higgs sector of the NMSSM consists of the Hu and Ha doublets of the MSSM along 
with a singlet S. The corresponding superpotential is given bjo 

WDXSHu-Hd+^nS^ (1.1) 

where A-B := e"''' AaBb with e^^ = +1. The soft supersymmetry breaking operators within 
the Higgs sector are 

-^soft D mljH^\^ + mljH,\^ + ml\S\'+(xA,H^-H,S + ^KA,S^ + h. A (1.2) 

The NMSSM Lagrangian as shown has a Z3 discrete symmetry that forbids the appearance of 
linear and quadratic singlet operators. When the singlet obtains a VEV in the early universe, 
this symmetry is broken spontaneously in one of three degenerate vacua, and dangerous 
domain walls can form [29]. These can be avoided by including a relevant operator that 
softly breaks the Z3 well below the electroweak scale. Such operators can arise in a natural 
way [301 EI] such that they eliminate domain walls, but are too small to have a significant 
effect on the Higgs boson phenomenology. 



^ We follow the notation and sign conventions of NMHDECAY 





S Hu Hd 


Q U^ D' L E^ 


PQ 

R 

PQ' 


-2 1 1 

2/3 2/3 2/3 
1 -2 


-10 0-10 

2/3 2/3 2/3 2/3 2/3 
0-12 0-1 



It is useful to introduce a complete basis for the set of all f/(l) transformations on S, 



;i.3) 



All of these couplings are broken by terms in the superpotential or the soft SUSY-breaking 
Lagrangian. U{1)pq is broken exphcitly by k and nA^, U{1)r is broken by XAx, kA,^, and 
the other trilinear soft terms and gaugino masses, while f/(l)pQ is broken explicitly by A 
and XAx. Throughout this paper, we will implicitly make a f/(l)'pQ field redefinition such 
that A is real and positive, as well as a U{1)pq transformation such that k, is real. With a 
U{l)ji rotation, we can take gaugino masses to all be real and positive provided they have no 
relative phases. The A terms subsequently generated by RG running will then be (mostly) 
real as well. Having fixed k to be real, but not necessarily positive, there is still a residual 
Z2 subgroup of U{l)pQ that can be used (in conjunction with f/(l)y), to make the VEVs 
(if°) = Vu and (H^) = Vd real and positive. The values of {S) = Vg and k can then take 
either sign. 



1.2 Minimization Conditions and a Challenge 

The Higgs-sector parameters in the NMSSM are {A, k, Ax, A^, fn'jj^, m?jj , wi|}. Electroweak 
symmetry breaking provides three minimization conditions that allow three of these param- 
eters to be exchanged for three VEVs. While it is a common procedure to solve for the soft 
masses m^^ and m^j in terms of v and tan /?, we find that it is instructive to instead solve 
for m|, /i = ^eff = Xvs, and k in terms of the other parameters. 

Demanding that f„, Vd, and Vg are all non-zero and real, we find 



171% — 171% tan^/9 1 

K ^ sin 2/? / AV rnl^+ml \ _ ^ 
A 2V/^2+ ^2 J 1^ 



2 + ^r+ 9 ' ]-—, (1-5) 



11^ A^ 2 /i A /i^ A /i 

These relations imply a difficulty for the NMSSM in gauge mediation [21] . We now review this 
argument. To leading order, the trilinear singlet soft terms and m| vanish at the messenger 
scale and are generated primarily from RG running to the electroweak scale. For lower 
mediation scales, the logarithm does not compensate the loop factor suppression, and these 
terms are expected to be somewhat smaller than the other soft terms. On the other hand 
under the condition that the lightest slepton is sufficiently heavy, the soft terms generated 
by minimal GMSB force \7n\j \ to be larger than about (250 GeV)^. Eq. (II. 4p then forces 
/i^ to be about the same size as \tTi\jJ\ (for moderate tan/3). These large values of /i^ lead to 



the hierarchies }? v^ j [x^ <^ 1 and \A\^f^l ii\ ^ 1. To satisfy Eq. (11.61) given the small values 
of m|//i^ expected from minimal gauge mediation, it is necessary to have |/t/A| ^1. As a 
result, sin 2/3 -C 1 is needed for Eq. (11. 5p to be satisfied. 

With both |k/A|, sin2/? -C 1, it is difficult to obtain an acceptable pattern of symmetry 
breaking. In this limit, the phenomeno logically interesting local extremum (corresponding 
to the solutions of Eqs. (I1.4fl.5lll.6l) where Vu-, Vd, and Vs are all non-zero) tends to be a 
saddle point rather than a minimum. This can be seen in the determinant of the CP-even 
Higgs boson mass-squared matrix. 

With |k/A| <^ 1, this expression is negative unless A^ <^ 1 is also tiny. Very small values of 
A and k, with /i fixed, correspond to a decoupling limit of the NMSSM in which the singlet 
states couple only very weakly to the rest of the theory, and the Higgs phenomenology reduces 
to that of the MSSM. The argument presented here applies to minimal gaugino mediation 
as well. 

Eqs. (Il.4til.6l) also indicate a way out of this difficulty: a deformation of minimal gauge 
or gaugino mediation that generates a large negative value of m| near the electroweak scale. 
In the present work, we investigate the phenomenology of two of the simplest possibilities for 
modifying the low-scale value of m^ in gauge and gaugino mediation. The first modification 
we consider consists of treating m| as a free parameter at the mediation scale of supersym- 
metry breaking. The second consists of adding charged vector-like states to the theory with 
superpotential couplings to the singlet ^7\ [221 121] • These couplings drive the singlet soft 
mass negative in the course of RG running, in a manner analogous to the effect of the top 
Yukawa coupling in the MSSM. The goal of the present work is to determine the viability 
of these modifications and to investigate the LHC Higgs signatures they predict. We do 
not dwell on a reduction of fine-tuning relative to the MSSM. Rather, our focus is to better 
understand under what conditions one might expect novel Higgs boson signatures. Then, if 
such signatures are indeed observed at colliders, we will have an important clue about the 
identity of the underlying theory. 

The plan of the remainder of this paper is as follows. In Section [2] we derive analytic 
expressions for the Higgs boson masses and couplings in the approximate ^(1)^ limit. In 
Section [3] we study the NMSSM Higgs sector within gauge mediation with an additional 
boundary contribution to m^ at the gauge messenger scale. Next, in Section H] we investigate 
the effect of adding charged vector-like states on the NMSSM Higgs states within minimal 
(unmodified) gauge mediation. We study analogous deformations of gaugino mediation 
within the NMSSM in Section [51 Section [61 is reserved for our conclusions. Some of our 
technical results are collected in Appendices \M and [Bl 



2 NMSSM Higgs Bosons with an Approximate U{1)r 

We begin by deriving approximate analytic expressions for the Higgs boson masses, mixing 
angles, and couplings in the NMSSM when the theory has an approximate U{1)r symmetry in 
the singlet sector p^ . This approximate symmetry is realized numerically as the hierarchies 
I^a/a^I -C 1 and Ivlre/Ail <^ 1- We shall also assume \Xv/fi\ <S 1, which is valid throughout 
most of the parameter space of the models we study in the coming sections. Further details 
about this expansion are listed in Appendix |Al 

To leading non-trivial order in the small ratios, and assuming k/A and sin 2/5 are not too 
terribly small, the tree-level CP-odd masses are 

2 / AAa\2A/K\2 i^ / „AAa\22 



The expression for m"^^ vanishes in the limit Ax^^ -^ showing that this state is the pNGB 
of the approximate U{1)r. This state is primarily singlet, while the A^ state is similar to the 
Higgs pseudoscalar in the MSSM. We are able to make this identification because the mixing 
among the singlet and the non-singlet pseudoscalars is suppressed by a factor of \v/fi <C 1. 
Full expressions for the mass and mixing matrices in this expansion are listed in Appendix 1X1 



Applying this expansion to the CP-even masses yields 

(2.3) 
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(2.4) 

(2.5) 

Among these states, h^ is SM-like, if ° is similar to the corresponding state in the MSSM, and 
hs is predominantly singlet. Again, this identification is possible because the mixing between 
the MSSM states and the singlet is suppressed by factors of Xv/fi <^ 1. The CP-even sector 
mixing matrices are also listed in Appendix Rl 

The (exact) mass of the charged Higgs is 

In the limit of /x^ ^ X^v"^, A\, A"^ this coincides closely with the masses of the A^ and H^ 

states. 

The coupling between the SM-like h^ state and pairs of the light a^ pseudoscalars is 
particularly important for the phenomenology of this scenario. This coupling corresponds 

6 



to the operator 

CD -^vh^aiai. (2.7) 

v2 

Within the expansion, the coefficient c is given by 

\A('--2.,,]('-^.,,)'^ (2.8) 



2 J \K J \^ /2/i 
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The ffist hne in this expression coincides with the resuh of Ref. |15], and corresponds to the 
axion-hke derivative couphng of the hght pNGB pseudoscalar. The second hne is new to our 
calculation. It is useful because it captures the contributions to the coupling that arise from 
the explicit breaking of the U{1)r symmetry by the A-terms. In terms of the coefficient c, 
the decay width for h^ — > astts is 

r(/i° -. a,a,) = -p^ (l - 4^) ' ' . (2.9) 

The approximate expression in Eq. (12. 8p agrees well with the full numerical result from 
NMHDECAY[28] in the appropriate limit. 



3 Gauge Mediation in the NMSSM with a Free m 



The first scenario we consider is a deformation of minimal gauge mediation in which the 
value of the singlet soft mass rrig taken to be a free parameter at the messenger scale M. All 
other soft terms are set to their standard gauge mediated values at scale M. Without this 
deformation, r7i|(M) vanishes at the leading order. This feature is the primary obstruction to 
merging gauge mediation with the NMSSM |21j. By liberating m'g from its minimal GMSB 
boundary condition we avoid this obstacle by fiat. The same deformation was considered in 
Ref. [27]. We expand upon and extend their results. 

We do not have a particular model in mind for how this deformation of minimal gauge 
mediation could arise. However, as a gauge singlet, one might imagine that S" is on a special 
footing and might feel the mediation of supersymmetry differently from the rest of the 
MSSMo A number of recent works have considered related modifications of gauge mediation 
within the NMSSM with this fact in mind [231 [25l EZ] • Prior explicit constructions typically 
generate new contributions to Ax and A^. This spoils the approximate U{1)r symmetry, 
and prevents the interesting decays to pseudoscalars. Modifying m| without altering Ax and 
A^ requires coupling S directly to a source that breaks supersymmetry but not the U{1)r. 



^This philosophy is similar in spirit to |321 133] . where supersymmetry boundary conditions for the Higgs 
fields were chosen to be different from the matter fields. 



We do not explicitly write a full model of supersymmetry breaking for the singlet field. 
However, we point out that supersymmetry breaking without spontaneous U{1)r breaking is 
a generic feature of simple O'Raifertaigh models [31], and that some care and complication 
is often required to break the U{1)r in this context [SU [35]. It is a model-building challenge 
to ensure that the additional SUSY breaking felt by the singlet field is of the same order as 
the rest of the SUSY breaking. 



3.1 Allowed Parameter Regions 

We begin by searching for phenomenologically acceptable regions of the NMSSM parameter 
space subject to GMSB boundary conditions for all soft terms other than ?n|. We specify 
the superpotential coupling A, as well as the supersymmetry breaking scale F/M, at the 
GMSB messenger scale, M. We specify tan/3 near the electroweak scale. With these inputs 
set, we compute the resulting low-energy spectrum. The soft supersymmetry breaking 
parameters at scale M, other than m|, are set to their gauge- mediated values assuming 
a minimal GMSB sector with a single set of 5 © 5 messengers. We use a modified version 
of NMSPEC/NMHDECAY [28] to perform the RG evolution of the model parameters, to 
find the low-scale values of k, \h\, and rri|, and to compute the low-energy spectrum and 
constraints. In general, the value of r7i| obtained in this way does not agree with the GMSB 
boundary condition of m|(M) ^ 0. Our scan encompasses the parameter ranges 

A G [0.001,0.7], tan/3 G [1,50], M G [10^ 10^^] GeV, F/M G [2,40] x 10^ GeV. (3.1) 

These ranges lead to low-scale gluino masses between about 350 GeV and 2500 GeV. 

The results of our parameter scans are shown in Fig. [T]within the A— ^ and A— tan f3 planes. 
The red (dark) points in these plots are consistent with all relevant phenomenological bounds 
except for the collider constraints on the neutral CP-even and CP-odd Higgs bosons, while 
the green (light) points satisfy the Higgs constraints as well. We do not demand that the 
lightest neutralino be the LSP since for gauge mediation the gravitino is usually expected 
to be the true LSP. However, we do require that the couplings remain perturbative up to 
the scale M, but not Mqut- This eliminates points with large A and k for larger mediation 
scales. In the limit M — > Mqut, we find that perturbativity requires (A^ + k^) < 0.45 near 
the electroweak scale [H [25]. When M < Mqut, additional charged states that enter the 
running at M, such as the messengers themselves, can help to slow down the growth of A 
and K above the messenger scale [36j . 

In the left panel of Fig. [H the lower limit on the red (dark) region corresponds to the 
condition det A1| > 0. For smaller values of A, this cutoff agrees with the relation given in 
Eq. (II. 7p . while more generally, it coincides with m^o > using the expression in Eq. (12. 3p . 
The upper boundary of the allowed red (dark) region in Fig. [1] can be understood from an 
examination of Eq. (11.51) . There is a close relationship between k/X and sin 2/3 when |/i| is 
much larger than A v and the singlet trilinear couplings. Since sin 2f3 is necessarily bounded 
by 1, so is k/A. 

It is clear from Fig. [1] that there are two disjoint populations of phenomenologically 
consistent (light/green) parameter points. The first population, which we call Region I, 




Figure 1: Allowed parameter points for the NMSSM with minimal GMSB boundary 
conditions for a single set of 5 © 5 messengers, but with no restrictions on m|. The red 
points do not include Higgs constraints while the green points do. 

has larger values of A > 0.4, k/A on the order of unity, and smaller values of tan/3 < 2.5 
(sin 2/9 > 0.7). The second population. Region II, has A < 0.08, k/A well below unity, and 
larger values of tan/5 > 5 (sin 2/3 < 0.38). For values of A between the large and small 
values taken on within the disjoint green regions, we find that the SM-like h^ Higgs boson 
is too light to satisfy the LEP bounds. In the larger A portion. Region I, the singlet F-term 
becomes important for smaller tan/5 and is responsible for increasing the mass of the h^. 
At smaller A, in Region II, 171^0 is very close to the value it would have in the MSSM, and 
larger values of tan/5 are necessary to push it above the LEP II bound. We discuss the 
phenomenology of Regions I and II below. 



3.2 Region I: Higgs Decays to Pseudoscalars 



Region I consists of points with A ~ k > 0.4 and tan/5 < 2.5. The Higgs phenomenology in 
this region can be very different from the SM when rUa^ < Tnho/2 and the /i° state decays 
predominantly into pairs of a^ pseudoscalars. Evidently this requires a light a^ pseudoscalar 
and a sizeable effective coupling c. When this is not the case, the Higgs phenomenology 
turns out to be very similar to the MSSM at large m^o. 

As discussed in Section [21 a light a^ pseudoscalar can emerge from the spontaneous 
breakdown of an approximate U{1)r symmetry in the singlet sector. This U{1)r is broken 
explicitly by the singlet trilinear A terms, but remains a good approximate symmetry 
provided they are much smaller than /i. In Fig. [2] we plot the values of Ax against A^, 
as well as Ax against fi for points within Region I. The green (light) points are consistent 
with all phenomenological bounds, while the blue (dark) points also have nia^ < m/io/2. 
From the right panel of this figure we see clearly a hierarchy between the singlet A terms 
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Figure 2: Values of Ax, A^, and n among phenomenologically acceptable parameter points 
in Region I described in the text. The green (light gray) points satisfy all coUider 
phenomenological constraints, while the blue (dark gray) points also have a light pseudoscalar 
with rria^ < m/,o/2. 



and /i. 

It is not hard to understand how the approximate U{1)r can give rise to rua^ < myip/2 
by comparing Fig. [2] with the results of Eq. (12. ip and Eq. (12.31) 1^1 To obtain such a light 
pseudoscalar, very small values of A^^ are necessary. There is a contribution to ma^ that 
goes as /i Afj, whereas the mass of the light h^ Higgs boson is largely controlled by \v in this 
region. Recall that }?v^ j [x^ -C 1. It is also necessary to have A\l [i reasonably small, although 
this requirement is much less severe due to the additional suppression by }?v^ j \x^ ^ 1 in 
Eq. dri]). 

Among the parameter points with a sufficiently light a^, the dominant contribution to 
the coupling c in Eq. (12.81) comes from the term involving A\. Thus, the small amount of 
[/(1)r breaking induced by the gaugino masses and transmitted to A\ in the course of RG 
running plays a dual role. It must be small enough to keep the pseudoscalars light, but still 
large enough to facilitate /i° — *> agag decays. As shown here, and previously observed in ^li] 
(there with breaking at the GUT scale), the size of the A terms derived from running can 
provide the right amount of U{\)r breaking to accomplish both of these tasks. 

The low-scale values of A\ and A^^ are generated in the course of RG running down 
from the messenger scale. They are sourced indirectly by the gaugino masses. Thus, the 
precise values of A\ and A^^ near the electroweak scale are sensitive to the GMSB parameters 
F jM and M. In Fig. [3] we show the values of A\ and A^ among allowed parameter points in 
Region I as functions of both FjM and M. We exhibit allowed points both with and without 



^ Numerically, we find the tree-level expression of Eq. (|2.ip to be accurate only for larger values of vria^ 
due to additional quantum corrections. However, Eq. (|2.ip remains useful for determining when the Og state 
will be light. 



10 



> 

o 



+ A. , all points 
A, , m < nt / 2 

. A . , all points 

. A . , m < m, II / 2 ** 













10 10 10 10 10 10 10 10 10 " 10 

M (GeV) 



1x10 2x10 3x10" 4x10" 

F/M (GeV) 



Figure 3: Dependence of the singlet trilinear Ax and A^ terms on tlie messenger scale M 
and tlie scale of supersymmetry breaking F/M. We exhibit phenomenologically allowed 
parameter points, as well as points satisfying the additional condition of rria^ < mho/2. 

rria^ < mho/2. Not surprisingly, smaller values of the trilinear A terms are obtained from 
lower values of F/M, corresponding to a lighter superpartner spectrum. For higher messenger 
scales M, the logarithmic enhancement of the singlet A terms becomes stronger, leading to 
larger values of these couplings. This tends to push up the mass of the a^ pseudoscalar, but 
it also helps to enhance the decay width of h^ —>■ agag allowing this mode to dominate over 
decays to bottom quarks. 

In Fig. m we show the masses of the lightest scalar h^ and the lightest pseudoscalar a^ 
for allowed parameter points in Region I in the left panel. In the right panel we show 
the dominant branching fractions for these states. We see that the decay properties of 
the h^ Higgs differ markedly from those of a SM Higgs when the branching fraction for 
h^ -^ agag is close to unity. When this is the case, the h^ state can be significantly lighter than 
114 GeV and still be consistent with the bounds from LEP II, as discussed in Ref. [10]. For 
pseudoscalar masses larger than about 10 GeV, the a^ decays primarily into bb, weakening 
the bound on the h^ mass to roughly 110 GeV. Very light pseudoscalars, below about 
10 GeV in mass, tend to decay mostly into rf (unless they are extremely light), allowing for 
h^ masses as low as about 90 GeV to be consistent with the bounds from LEP II pjj. In 
the present context, it is challenging to make A^ small enough to get such a light a^. We 
find that rua^ < 10 GeV generally requires a cancellation between the A^ and Ax terms in 
Eq. (12. ip . implying a degree of fine-tuning. This possibility is constrained by and will be 
further probed by searches for T —^ 0^7 [371 EH] as well as measurements of {g — 2)^ [39]. 

While the decays of the SM-like h^ Higgs boson can be significantly modified by the 
presence of a light singlet pseudoscalar, the other properties of the h^ state such as its 
production rates and decay widths into SM final states generally remain nearly unchanged 
relative to the SM Higgs. For example, the mixing between the SM-like combination of the 
CP-even gauge eigenstates {i.e. the combination that has the same couplings with the gauge 
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Figure 4: Masses and branching fractions of the hghter h^ and a^ Higgs bosons for allowed 
parameter points in Region I described in the text. 



bosons as in the SM) with the singlet in the h^ is given by the matrix element 

XXvfX 

Ul2 ^ S2/3 

2k jj, \K 



(3.2) 



Here we have kept only the leading terms in the expansion in A v/fi and Ax.^./fJ'- More details 
about the mixing matrix are given in Appendix |Xl Numerically, we find |?7i2| ^ 0.10, along 
with \Uii\ > 0.995. In particular, the dominant Higgs production channels via loop-induced 
coupling to gluons or direct couplings to gauge bosons will only be reduced by a factor of 
Ufi, or less than about 1%. 

The LHC signatures of the h^ Higgs depend strongly on its branching fraction into pairs 
of Gg. If these decays are suppressed, either by kinematics or a small coupling, the Higgs 
signatures will be very similar to the MSSM. When BR{h^ -^ OsOs) is non-trivial, the 
branching fractions of the /i° into MSSM final states will be reduced according to 



BRi = BRf'^'^^^ [1 - BRih'^ 



aM^ 



Despite this suppression, some of the more promising SM Higgs channels such as h^ 



(3.3) 



77 



[mho 



< 



125 GeV) and /i° -^ WW* (m^o > 125 GeV) could still be visible in their own right 



even with BR{hP 



a.M. 



close to its maximal value of about 0.9. 



It may also be possible to search for the h^ Higgs boson through its decays to light 
pseudoscalars if both BR{hP — >• astts) and BR{as — ^ bb) are close to unity [IHIHT]. This search 
relies on Higgs production in association with a weak gauge boson to reduce the background. 
It is challenging because it requires multiple 6-tags to reduce the background, and therefore 
requires a high 6-tag efficiency along with a good understanding of the mis-tagging rate. 
Even so, Ref. [4T| finds that it should be possible to discover the h^ Higgs through this 
channel with at least 30 fb~^ of integrated luminosity. When the a^ pseudoscalar is so light 
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that it decays predominantly into rf^ which we find to be fairly unlikely in Region I, Ref. [12] 
proposes a LHC search strategy using forward proton tagging (central exclusive production) 
with the proposed FP420 detectors |l3j. The recent study of Ref. [H] suggests that this 
decay channel might also be observable at the LHC in the weak boson fusion production 
channel. 

We have concentrated so far on the phenomenology of the h^ Higgs boson and the a^ 
pseudoscalar. The other Higgs states, the CP-even H^ and hs along with the CP-odd AP and 
the charged if^, could also turn up at the LHC. However, among the phenomenologically 
allowed points in Region I these states are all heavier than about 500 GeV and mix only very 
weakly with the SM-like h*^ . This makes them difficult to detect at the LHC. Since tan/? is 
small in Region I, these states have significant branching fractions for decays into ti. The hg 
scalar can also decay efficiently into pairs of a^ pseudoscalars, while the A^ pseudoscalar can 
also decay appreciably to the lighter superpartners. Given the large masses of these states 
along with their weak production cross-sections, we do not expect that they will produce a 
significant signal at the LHC [l5l B6]. The charged H^ Higgs decays primarily into t6, as 
well as h'^W^ , and may possibly be visible through the former mode j45l HGj. 

Besides the new Higgs bosons, the singlet superfield S also gives rise to an additional fifth 
neutralino state. For the parameter values in Region I, this state consists primarily of the 
fermion component of 5* with a small higgsino admixture, and has a mass close to 2k^/\. As 
k/A ~ 1 and \^\ > Mi 2 in Region I, this state is typically the heaviest neutralino. Moreover, 
its small mixing with the MSSM neutralinos implies that it will have a low production rate 
at the LHC, either by direct creation or through SUSY cascade decays. Cascade decays 
could also provide an interesting source of h^ Higgs bosons, although we do not pursue this 
possibility here. 



3.3 Region II: (Partial) Singlet Decoupling 

This region corresponds to the phenomenologically allowed points in Fig. [1] with A < 0.1, 
k/X < 0.4, and tan/5 > 5 (sin 2/5 < 0.4). As in Region I, the minimal GMSB boundary 
conditions lead to values of |/i| that are much larger than Xv and the singlet trilinear Aa,^ 
couplings. This signals a certain amount of fine-tuning in the electroweak symmetry breaking 
conditions. From the above values and the minimization conditions of Eqs. (IL4IIL5|1.6p . we 
see that this region has |m|//i^| <^ 1 near the electroweak scale. The small values of A ^ 1 
that arise in Region II are needed to satisfy the LEP bound on the SM-like Higgs boson 
mass. Expanding Eq. ( 12. 3p in small n/X and sin 2/5, we find that 

AV(<l^-^V ,3,4) 



,2 _ \2„,2 



2A2 K 



To obtain a sufficiently large Higgs mass, it is then necessary that (?^/2A^ ^ }? j k^ ^ forcing 
A2 « 1. 

The Higgs sector in Region II always contains a very light pseudoscalar a^ as well as a SM- 
like /i° Higgs boson. The pseudoscalar mass is protected by an approximate U{X)r symmetry 
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Figure 5: Values of Ax, A^, and fi among phenomenologically acceptable parameter points 
in Region II described in the text. 



in the singlet sector due to the GMSB boundary conditions, in addition to an approximate 
U{l)pQ arising from the small value of n/X. On account of this double-protection, the 
electroweak scale value of A^ is particularly small among the allowed parameter points. 
This is clearly shown in the left panel of Fig. O From the right panel of Fig. [5] we also see 
that the value of \fi\ is much larger than A^, A\, and Xv within the allowed parameter space. 

With these very small values of Ax and A^,, the mass of the pseudoscalar a^ in Region II 
is always much less than half the mass of the SM-like h^ Higgs boson. This feature is 
illustrated in the left panel of Fig. [61 In the right panel of Fig. [6] we show the dominant 
branching fractions of the h^ Higgs boson. Despite being kinematically allowed, we find that 



h' 



a^a. 



almost never occurs, with BR{h^ 



a.a. 



< 3 X 10 . Instead, the decay fractions 



of the h^ Higgs are nearly identical to those of a SM Higgs. The reason for this suppression 



oih^ 



aM. 



can be seen in Eq. (12.81) : the coupling c is proportional to A^ -C 1. 



Collider production of the h^ Higgs in Region II is also nearly identical to the SM. The 
singlet component of this mass eigenstate is quite small, corresponding to the mixing element 
|f^i2| ^ 0.25, while the suppression of Higgs production through the usual channels is down 
by only Ufi > 0.94. (See Appendix |A] for a full account of the mixing matrices.) Thus, the 
LHC signatures of the h^ in Region II will be nearly identical to those of a SM Higgs (or an 
MSSM Higgs at large M^o). The mixing among the CP-odd states is even more suppressed, 
with the light as being almost exclusively singlet. This provides another way to understand 
why the branching fraction of the h^ into pseudoscalar pairs is so small in Region II. 

Among the heavier Higgs bosons, the mostly- singlet hg can be relatively light since its 
mass is suppressed by a factor of k/\ <C 1. We find masses as low as 170 GeV, and as high 
as 1500 GeV. While this state is mostly singlet, a mixing with the SM-like Higgs as large 
as If/21 1 — 0.25 is possible. The primary decay modes are WW, ZZ, and h'^h^. For larger 



mixings and moderate masses, in the range 200 GeV 



< 



mhs 



< 



500 GeV, the hg may be 
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Figure 6: Masses of the lighter /i° and a^ Higgs bosons, and the dominant branching fractions 
of the h^ Higgs for allowed parameter points in Region II described in the text. 

visible at the LHC through its ZZ final state. All the other Higgs bosons, the H^, A^, and 
H^, have masses in excess of 800 GeV, and are challenging to find at the LHC. 

The rest of the particle spectrum within Region II is similar to the MSSM, but with 
an additional neutralino from the fermion component of S. To push ruho above the LEP- 
II bound, the superpartner spectrum must be somewhat heavy with M3 > 800 GeV. On 
the other hand, the corresponding mostly-singlet neutralino state state can be as light as 
170 GeV in Region II, and can even be the lightest superpartner aside from the gravitino. 
Due to the very small values of A and k, this state is almost pure singlet with a tiny 
higgsino component, and couples only very weakly to the rest of the MSSM. It will therefore 
almost completely decouple from the LHC phenomenology unless it is the NLSP (with a 
gravitino LSP). In this case, the presence of an additional state can modify cascade decay 
chains [171 HHl 139] and possibly also give rise to neutral displaced vertices from decays of the 
NNLSP to the NLSP ^47,i50j. A mostly singlet NLSP could also be problematic if it decays 
after the onset of nucleosynthesis [511 |52l [53] . While a full analysis of this issue is beyond 
the scope of the present work, we expect that lighter gravitino masses (leading to shorter 
decay times) might be necessary to ensure that these decays occur before the light elements 
are formed. 



4 Gauge Mediation in the NMSSM with Exotics 

As a second variation on minimal gauge mediation in the NMSSM, we consider adding 
charged vector-hke exotics to the theory [171 [211 [22] . We study exotics in the form oi DQ)D'^, 
and L © L"^, with D = (3, 1, —1/3) and L = (1,2, —1/2). Taken together, these exotics have 
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c 



the quantum numbers of a 5 © 5 of SU{5). We include the trihnear superpotential couphngs 

Wd^dSDD' + ^lSLL'', (4.1) 

as well as a corresponding set of soft supersymmetry breaking operators. These couplings 
tend to drive m'g negative in the course of RG running from the gauge messenger scale 
down to near the electroweak scale, thereby facilitating singlet condensation and electroweak 
symmetry breaking. The relevant RG equations for this evolution with Nd sets oi D (B D 
and Nl sets of L © L"^ are listed in Appendix [Bl Note that we assume the exotics do not act 
as gauge messengers, which can be enforced with the approximate Z3 symmetry discussed 
in the introduction. 



4.1 Allowed Parameter Regions 

We have searched for regions of the NMSSM parameter space that lead to an acceptable 
phenomenology by scanning over the model parameters. Our strategy consists of specifying 
^j ^D, ^L, and F/M at the GMSB messenger scale M, along with tan/5 near the electroweak 
scale, and computing the low-energy spectrum that results from these input parameters. 
In doing so, we use a modified version of NMSPEC/NMHDECAY [28j to perform the RG 
evolution of the model parameters, to find the low-scale values of k, \h\, and m|, as in 
Eqns. (11.41 II. 5111. 6p . and to determine the phenomenological constraints. The value of ?ti| 
computed in this way will not usually agree with the GMSB boundary condition of r7i|(M) ~ 
0. To correct for this, we adjust the value of ^^ (or C,l), repeat the running, and iterate until 
the value of |m|| at the input scale M lies below a small cutoff value. The parameter ranges 
covered in our scans are 

A €[0.001,0.7], eLe[o,i], eDe[o,i], 

(4.2) 
tan/? e [1, 50], Me [10^ lO^^] GeV, F/M e [2, 40] x 10^ GeV. 

We assume minimal GMSB boundary conditions for all the soft terms with a single set of 
5 © 5 messengers. 

In Fig. [7] we show the allowed parameter points obtained by scanning with a single set 
A/'s = 1 of 5 © 5 (non-messenger) exotics. In the left panel of this figure we show points in 
the X-K plane, while in the right panel we exhibit points in the A-tan/5 plane. The red (dark) 
points agree with all relevant phenomenological bounds except for the mass constraints on 
the neutral CP-even and CP-odd Higgs bosons, while the green (light) points are also 
consistent with the Higgs constraints. As in the previous section, we do not demand the 
lightest neutralino be the LSP: in gauge mediation the gravitino is usually expected to be 
the true LSP. We do require that the couplings remain perturbatively small up to the scale 
M (but not Mgut)- 

The phenomenologically acceptable parameter points in Fig. [7] are similar to the pa- 
rameter points found in Region II discussed in the previous section. These points all have 
very small values of A, and relatively small values of k/A and sin 2/3. Small values of A in 
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Figure 7: Allowed parameter points for the NMSSM with minimal GMSB boundary 
conditions for a single set of 5 © 5 messengers, and a single set of 5 © 5 exotics. The 
red (dark) points do not include Higgs constraints while the green (light) points do. 

this region are needed to make the SM-like Higgs boson sufficiently heavy. Unlike in the 
previous section, however, we do not find any points at larger values of A and n, analogous 
to Region I . For points of this type to lead to an acceptable pattern of symmetry breaking, 
Eqs. (I1.4|1.5|1.6p require |m||//i^ ~ 1 (or large A terms). While the new Yukawa couplings 
in Eq. (14.11) help to drive m| to more negative values, their effect is not strong enough to 
open a region of parameter space with A ~ k ~ 1. With a single set of non- universal exotics, 
such as a lone L (B L'^ or D (B -D^, we find qualitatively similar results. 

To magnify the effect of the exotics on the running of m|, we have also performed scans 
with multiple sets of 5 © 5 (non-messenger) exotics. For simplicity, we assume universal 
values of the couplings ^d and ^l for all flavors of exotics. Adding even only a second set of 
5 © 5 exotics opens a new region of phenomenologically consistent points, similar to Region I 
discussed in Section |3l These points have «: ~ A > 0.4 and tan/3 < 2.5. A similar region 
can emerge with more than two sets of 5 © 5 exotics. The allowed region with A <^ 1 also 
remains. With several sets of exotics, a tension arises between larger values of A and k at the 
low scale and perturbativity up to the messenger scale M since the exotic Yukawa couplings 
of Eq. (14. ip speed up to the running of A. 

The NMSSM Higgs boson spectrum depends sensitively on the values of Ax, A,^, and fi 
near the electroweak scale. In Fig. [8] we plot these parameters for A^s = 1 and N^ = 2 sets 
of 5 © 5 exotics. For the case of A^s = 2, we split up the points according to whether they 
fall into the A < 0.1 region, or the A > 0.4 region. These plots indicate that A\ and A^, 
remain smaller than |/i| (larger values of Aa,^ correspond to very large values of /x), but are 
enhanced relative to the electroweak scale. This enhancement comes from the contributions 
from the new Yukawa couplings ^d and ^l in Eq. (14. ip to A\ and A^ in the course of RG 
running. 
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Figure 8: Values of Ax, A^, and /i among phenomenologically acceptable parameter points. 
We exhibit allowed points with N^ = 1 and A < 0.1, as well as N^ = 1 and A < 0.1, and 
A^5 = 2 and A > 0.4. 



4.2 LHC Higgs Signatures with Exotics 

Much like in Section [3l the Higgs phenomenology with additional vector-like exotics can be 
split into two classes, depending on whether A is small (A < 0.1) or A is larger (A > 0.4). 
This second possibility only occurs with at least two sets of 5 © 5 exotics. We will discuss 
both of these cases in turn. 

In the allowed parameter region with A ^ 1 (along with \k/X\ ^ 1 and sin 2/5 <^ 1) 
the Higgs boson phenomenology is quite similar to what we found in Region II discussed 
in Section [3l We show the masses of the neutral CP-even and CP-odd Higgs bosons in 
Fig. [9] for phenomenologically consistent points obtained in our scan with a single set of 
5 © 5 exotics. The lighter Higgses consist of the SM-like /i°, the mostly-singlet hs, and a 
light mostly-singlet pseudoscalar a^. The A^, H^, and H^ states are similar to their MSSM 
counterparts, and are generally very heavy. There are some slight but important differences 
compared to Region II of Section [3l however. 

The singlet A terms that arise in the course of RG running are now somewhat larger, due 
to the new exotic Yukawa couplings (see Fig.[H]). The result is the mostly singlet pseudoscalar 
as can be considerably heavier than in Region II discussed above. Depending on the values of 
these trilinear couplings, m^^ can range from below 1 GeV all the way up to over 250 GeV. 
This state is nearly pure singlet with only a tiny admixture of the A^ as a result of the very 
small values of A and n. We expect it to be almost completely decoupled from the rest of the 
theory, and consequently invisible at the LHC. Similar considerations minimize the ability 
to probe this state through upsilon decays, if kinematically accessible. The dominant decay 
modes of the a^ are into bb, tt, and so on, depending on which channels are kinematically 
accessible. 



1000 - _^ ***^^^^^^^^^:.\ 



o ; <;f • • 
s - -i' 

S -.:* , 

10- h 

^ h" 

: • H° 





' 1 


-I 


**+""■ 




- 


■ ^ 




-I 


i*' ' 




- 


: .-v ■ ' 




a 


- 


- ,'" 




a" 


- 








.1.1.1.1 


1 





> 

100 u 

o 



0.05 0.1 0.15 0.2 0.25 0.3 0.05 0.1 0.15 0.2 0.25 0.3 

Kl X kI X 



Figure 9: Neutral Higgs boson masses for phenomenologically allowed points in the NMSSM 
with a minimal GMSB spectrum and a single set of 5 © 5 exotics. 



The LHC phenomenology of the h^ and hs states is more interesting. As in Region II of 
the previous section, the branching fractions of the h^ Higgs are similar to those of a SM 
Higgs, and thus the LEP bound of about 114 GeV applies throughout most of the parameter 
region. An exception occurs when the mixing between the h^ and hg states is enhanced when 
they become nearly degenerate. Even in this case, the reduction in the mass bound on the 
SM-like state is usually less than a few GeV. The mostly-singlet hs Higgs can be very light 
in some cases, with masses as low as 20 GeV, but also as heavy as about 800 GeV. This 
state inherits its couplings to the MSSM through its mixing with the SM-like Higgs. Thus, 
the hs decay modes mirror those of a SM Higgs boson. The chief exception to this occurs 
when the channel hs — > h^h^ opens up. The corresponding branching fraction can be as large 
as BR{hs -^ h^h^) ~ 0.4. For heavier hg masses, this state might be visible at the LHC 
through its decays to ZZ, even though its production cross-section is suppressed relative to 
a SM Higgs by its large singlet component. All the other Higgs boson states are heavier than 
about 500 GeV, making them challenging to find at the LHC unless tan/? is very large. 

The rest of the particle spectrum within this small A region is similar to that of the 
MSSM. Compared to Region II, the overall superpartner scale can be somewhat lower, with 
gluino masses as small as M3 ~ 600 GeV now possible, and the additional mostly-singlet 
neutralino can be very light. The mass of this singlet neutralino is approximately 2Kfi/X, 
which can be as light as about 20 GeV for particularly small values of k/A, but also as large as 
several hundred GeV in other portions of the allowed parameter space. Mixing between the 
singlet fermion and the other neutralinos is always very small. The dominant contribution 
comes from the higgsinos with a mixing angle jOisI ^ 0.01. Thus, this fifth neutralino state 
is nearly invisible at colliders unless it is the NLSP (with a gravitino LSP). As discussed in 
Section [3], a mostly-singlet NLSP can modify sparticle cascade chains [171 HH |19]. In this 
case, lighter gravitino masses may be required to ensure that the NLSP decays safely before 
nucleosynthesis. 
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The second allowed region of NMSSM parameter space with heavy exotics and minimal 
GMSB boundary conditions has A > 0.4 and smaller values of tan/? < 2.5. These parameter 
ranges are similar to those encountered in Region I of Section [31 Despite the similarity of 
these values, the Higgs boson phenomenology does not mirror that of Region I. In fact, it 
does not significantly differ from that of the MSSM. To populate this large A region, at least 
two sets of 5 © 5 exotics are required in order to drive m^ sufficiently negative over the 
course of the RG evolution. The exotic Yukawa couplings, ^d and ^l, responsible for doing 
so also contribute to the low-scale values of Ax and A^. This ruins the approximate U{l)Ft 
symmetry in the singlet sector. As a result, the singlet pseudoscalar is no longer very light, 
with masses above ma^ ^ 350 GeV. The SM-like /i° state is therefore very SM-like, both 
in its production and decay modes, since it is no longer able to decay into pairs of the a^ 
pseudoscalar. All the other Higgs boson states are heavier than about 500 GeV, and will 
be difficult to find at the LHC, particularly with the smaller values of tan/? that occur in 
this region. The additional neutralino state that arises from the fermion component of the 
singlet S is also quite heavy, with a mass of 2 k/x/A ~ 2//, and has only a small mixing with 
the higgsinos. It too will be essentially invisible at the LHC. 



4.3 Phenomenology of the Exotics 

The vector-like exotic states we have added to the theory to facilitate electroweak symmetry 
breaking can themselves be a source of new signatures at the LHC. These charged exotics 
must certainly be heavy enough to have avoided detection already. Beyond this, the exotics 
are problematic for cosmology if they are overly long-lived (or stable). We briefly consider 
the additional bounds and potential signatures that arise from the exotics. These should 
inform any model building attempt. 

In Fig. [in] we plot the masses of the charged D- and L-type exotics originating from 
the superpotential couplings of Eq. (14.11) . As before, the red (dark) points in this plot are 
consistent with all phenomenological bounds other than the constraints on the Higgs sector 
(and on the exotics themselves), while the green (light) points satisfy all relevant Higgs 
bounds as well. Nearly all the points in this flgure have exotic masses well in excess of the 
current limits. 

The precise collider bounds on the exotics depend on whether they are long- or short- 
lived on collider timescales. In the case of long-lived charged leptons, the best bound comes 
from searches by OPAL for (effectively) stable charged particles, and is rrii > 100 GeV for 
a heavy spin-1/2 particle with electric charge ±1 and no color [M]- The bound for decaying 
charged leptons is about the same. The most stringent bound on long-lived heavy quarks 
comes from Tevatron searches for charged massive particles. A preliminary analysis with 
1 fb~^ of Tevatron Run II data suggests a limit on the cross section for such states of about 
0.1 pb [55j, corresponding to quark masses up to nearly 300 GeV. For short-lived heavy 
quarks, the precise bound depends on how they decay. CDF has performed a preliminary 
search for a heavy exotic top quark decaying through t' —>■ Wq, and flnd mf > 258 GeV [56] . 
No 6-tag is used in the analysis so we expect this result to apply for a short-lived exotic D 
as well [57] . 
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Figure 10: Masses of the charged D and L exotics among the phenonieno logically allowed 
parameter points found for the NMSSM with a single set of 5 © 5 exotics and a minimal 
GMSB soft mass spectrum. 

To be cosmologically acceptable, the heavy exotics cannot be too long-lived. If stable on 
the lifetime of the universe, heavy charged exotics are very stringently constrained by searches 
for anomalously heavy atoms. These bounds are so severe that even the tiny density of heavy 
exotics created by cosmic rays is unacceptably large [5B]. Long-lived heavy exotics are also 
dangerous if they decay after nucleosynthesis as their decay products can modify the light 
element abundances, distort the CMB blackbody spectrum, or contribute to cosmic rays [SH]. 
On the other hand, a significant coupling between the exotics and the MSSM matter fields 
can give rise to too much flavor mixing. 

In the NMSSM with an approximate Z3 discrete symmetry, the exotic D^^^ and L^^^ states 
can decay through the d = 4 superpotential operators 



LL'H^Hd, LL'SS, D'DHuHa, D'DSS. 



(4.3) 



Such operators can be consistent with the approximate Z3 symmetry while still allowing 
all the standard NMSSM operators, as well as the neutrino mass operator (LHu)^, and 
forbidding mixing between the exotics and the MSSM matter states at the renormalizable 
level. Decays through the operators of Eq. (14.31) occur safely before nucleosynthesis pro- 
vided the heavy mass scale suppressing them is less than about the GUT scale, Mqut — 
IQ^^ GeV [591 EOl. 



If the charged and colored exotics are not too heavy, they might lead to observable 
signatures at the LHC. Stable (on collider times) heavy quarks were studied recently in 
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Ref. [59]. These can form charged exotic hadrons that punch through to the muon chamber. 
By measuring the time of flight, they can be distinguished from ordinary muons. A significant 
signal of ten events with almost no background will be generated with 10 fb~^ for heavy 
quark masses below tjid < 1700 GeV, and m^ < 1450 GeV for the scalar superpartner. 
The precise limits on unstable heavy quarks depend on how they decay, but are generally 
of the similar size. Heavy stable charged leptons can also be detected at the LHC through 
time-of-fiight techniques [611 [62]. Masses below about 950 GeV can be detected in this way 
with 100 fb~^ of integrated luminosity. 



5 Gaugino Mediation in the NMSSM 

Gaugino mediation of supersymmetry breaking shares many of the attractive features of 
gauge-mediated supersymmetry breaking [191 EQ] . In the context of obtaining a light pseu- 
doscalar in the NMSSM, minimal gaugino mediation has the helpful property of vanishing 
trilinear A terms in the singlet sector at the mediation scale. Unfortunately, like minimal 
gauge mediation, this mediation mechanism in its minimal form has trouble generating a 
sufficiently negative low-scale value of the singlet soft mass r7i| to obtain an acceptable 
pattern of electroweak symmetry breaking. In this section we study gaugino mediation in 
the NMSSM with the same modifications of the previous two sections: we augment the 
theory by additional contributions to the input value of 7tt,|, or by adding new vector-like 
exotics coupled to the S field. 

5.1 Higgs Phenomenology with m| Free 

We study first minimal gaugino mediation within the NMSSM with the singlet soft mass 
m| taken to be a free parameter. One might motivate this scenario by putting the S field 
in the bulk, so that it could feel the SUSY breaking directly. However, to avoid generating 
singlet A terms that would ruin the approximate U{1)r symmetry in the singlet sector, such 
a bulk singlet should couple to a SUSY breaking source that also preserves the U{1)r. It 
is a model-building challenge to obtain such a source of supersymmetry breaking with the 
correct mass scale. 

As in Section [3], we search for phenomenologically acceptable parameter regions by 
scanning over the model parameters. Our scan encompasses the ranges 

A e [0.001,0.7], tan/3 G [1,50] GeV, (5.1) 

Ml/2 e [100, 1500] GeV, M^ G [10^ 10^^] GeV. 

The parameter M1/2 defines the gaugino mass at the compactification scale Mc through the 
relation 

M,{M,) = 2gl{M,)My2. (5.2) 

This choice corresponds to universal gaugino masses equal to M1/2 when M^ = Mqut 
{gl^McuT) — 1/2). The input values of all the soft scalar masses (save m|) and all the 
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A terms are taken to vanish at the input scale Mc- Over the course of the RG running, we 
require that all couplings remain perturbatively small up to M^. 

The allowed regions found in our scans for minimal gaugino mediation with rn^ free are 
nearly identical to what we found in Section [3] for minimal GMSB, and illustrated in Fig. [71 
As before, there are two disjoint phenomenologically consistent regions: one at small values 
of A, k/X, and sin 2/3 much like in Region II; and the other with larger values of A and 
sin 2/? ~ k/X ~ 1 similar to Region I. In both allowed regions, there exists an approximate 
U{l)ji symmetry in the singlet sector that is realized as a hierarchy between Ax and A,^ 
relative to /i. The Higgs boson phenomenology within these two regions is qualitatively the 
same as in Regions I and II discussed in Section [31 

Let us also mention that in the parameter scans, we do not demand that the lightest 
NMSSM superpartner be neutral. For compactification scales below about M^ < 10^^ GeV, 
the lightest MSSM superpartner in minimal gaugino mediation is often a mostly-right-handed 
stau [m [201 E3]. The true LSP in this case is typically the gravitino [HI [M]- At larger 
values of the compactification scale, the lightest MSSM superpartner in minimal gaugino 
mediation is usually a mostly-Bino neutralino. Within the NMSSM, there can also arise a 
very light mostly- singlet neutralino state when k/X is small, as we discussed in Sections [31 
and m Such a state can supplant the stau as the lightest NMSSM superpartner in gaugino 
mediation, providing another way to get around the problem of a charged LSP. This also 
leads to new possibilities for the production and identity of the dark matter, such as by 
the decoupling of a singlet neutralino LSP [19], or through the superWIMP scenario [51] 
with a singlet neutralino NLSP. This latter possibility, however, is likely constrained by the 
overproduction of hadronic debris after the onset of nucleosynthesis [5T| [53] . 



5.2 Higgs Phenomenology with Vector-Like Exotics 

A second modification of minimal gaugino mediation that can improve the prospects for 
electroweak symmetry breaking in the NMSSM is to add vector-like exotics to the theory. 
As in Section [H we consider 5 ©5 exotics with superpotential couplings to the singlet given by 
Eq. (14.11) . We again search for phenomenologically acceptable parameter regions by scanning 
over the model parameters. Our scans encompass the same ranges of input values as listed 
above, along with the exotic couplings 

^^€[0,1], ae[o,i]. (5.3) 

When we include more than one set of exotics, we assume that the values of C,d and C,l are 
the same for all exotic flavors. Gaugino mediated boundary conditions are imposed at the 
compactification scale M^ for all the soft terms, including m|. 

Our search strategy is similar to Section HI and consists of specifying A, ^d, ^l, and 
Ml/2 at the compactification scale M^, along with tan/3 near the electroweak scale, and 
computing the low-energy spectrum that results from these input parameters. In doing so, 
we use a modified version of NMSPEC/NMHDECAY [2ij to perform the RG evolution of 
the model parameters and to find the low-scale values of k, \fi\, and m|. The value of m| 
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computed in this way will not usually agree with the gaugino mediation boundary condition 
of m|(Mc) ~ 0. To correct for this, we adjust the value of ^^ (or ^l), repeat the running, 
and iterate until the value of |m|| at the input scale M^ lies below a specified small cutoff 
value. Again, we demand that all couplings remain perturbative up to Mc, but we do not 
require the lightest NMSSM superpartner to be neutral. 

The results of our scans are very similar to what we found Section IH and illustrated 
in Fig. [71 With a single set of 5 © 5 exotics, the phenomenologically allowed region of the 
parameter space is nearly identical to the small A region discussed in Section H] as well as 
Region II studied in Section [31 With two or more sets of 5 © 5 exotics, we find a second 
disjoint allowed parameter region with A ~ k > 0.4 and tan/3 < 2.5. This region of the 
parameter space is much the same as the large A region discussed in Sectional In particular, 
the new exotic Yukawa couplings C,d and C,l help to transmit the ^(1)_r breaking from the 
gauginos to the singlet sector, leading to somewhat larger values for the trilinear A terms 
at the low scale. The mostly-singlet pseudoscalar a^ is always heavier than the h^ Higgs as 
a result, leading to a Higgs collider phenomenology nearly identical to the MSSM with a 
pseudoscalar mass in excess of 500 GeV. 



6 Conclusions 

In the present work we have studied two simple deformations of gauge and gaugino mediation 
within the NMSSM. The first deformation consists of allowing the singlet soft mass m| to 
be a free parameter at the mediation scale, such as might arise if the singlet couples to 
a U{1)r preserving source of supersymmetry breaking. The second deformation involves 
adding vector-like exotics with superpotential couplings to the singlet superfield. Both 
deformations facilitate electroweak symmetry breaking by driving the singlet soft mass m| 
to negative values in the infrared. 

Near the electroweak scale, for either deformation of both minimal gauge and gaugino 
mediation, we find a hierarchy between the value of the effective /i parameter, /i = fieff = 
Xvs, and the Ax and A^ soft trilinear couplings, as well as Aw. The relative smallness of 
the singlet A terms leads to an approximate t/(l)i? symmetry in the singlet sector, and 
consequently to a light pNGB pseudoscalar when this would-be symmetry is spontaneously 
broken. The presence of this light pseudoscalar can have a large effect on the Higgs signatures 
of the theory. 

We find two distinct ways in which these deformations can lead to consistent electroweak 
symmetry with a phenomenologically acceptable Higgs boson spectrum. The first and more 
interesting case requires |m|//i^ | as well as k/X and sin 2/5 to be on the order of unity near the 
electroweak scale, along with A > 0.4. With these parameter values, there is an additional 
F-term contribution to the mass of the lightest CP-even Higgs boson. If, in addition, A^ 
and Ax remain small near the electroweak scale, there is a SM-like Higgs boson state in the 
spectrum that can decay predominantly into pairs of the light pNGB a^ pseudoscalar leading 
to new Higgs boson signatures at the LHC. These pseudoscalars usually decay primarily into 
bb, but can have a dominant branching into tt if they are particularly light {ma^ < 10 GeV). 
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This scenario is realized most easily in both gauge and gaugino mediation when ?7i| is allowed 
to be free, with no additional contributions to the singlet trilinear couplings at the mediation 
scale. 

The second way to obtain consistent electroweak symmetry breaking within the defor- 
mations considered here is to have k/A ~ sin 2/3 -C 1 near the electroweak scale. This occurs 
when |m||//i^ -C 1 with small singlet A terms. For the electroweak symmetry breaking 
extremum to be stable, it is then necessary to have A -C 1. In this case, the singlet sector 
couples only very weakly to the MSSM states, and therefore mostly decouples. It is still 
possible to have a very light pNGB pseudoscalar, but since it interacts only feebly with the 
SM-like Higgs state in the spectrum, the decays of this Higgs boson into pseudoscalar pairs 
are extremely rare. In general, the Higgs phenomenology in this case is very similar to the 
MSSM with a relatively heavy pseudoscalar Higgs boson. Despite the decoupling property of 
this region, the supersymmetric phenomenology can still be modified if the lightest NMSSM 
superpartner state is a mostly-singlet neutralino. 

Let us emphasize that we have only considered simple deformations of minimal gauge 
and gaugino mediation. In these minimal versions, we always find relatively large values for 
the effective /i parameter near the electroweak scale. By allowing for non-minimal versions of 
these mediation mechanisms, it may be possible for much smaller values of /i to emerge. In a 
scenario with /x ~ A f , there may be new ways to obtain a consistent pattern of electroweak 
symmetry breaking, possibly also with decays of a SM-like Higgs bosons into pairs of light 
pseudoscalars. While not all models of supersymmetry breaking will give rise to modified 
Higgs boson phenomenology, it is exciting that very minor modifications to models as simple 
and well-motivated as gauge mediation and gaugino mediation can. This emphasiszes the 
need to search for Higgs bosons at the LHC with as broad a net as possible. 
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Appendix 

A Higgs Boson Masses and Mixings 

In this Appendix we present analytic expressions for the Higgs boson mass eigenvalues and 
mixing matrices valid in the limit of an approximate t/(l)_R symmetry in the singlet sector 
and Xv <ti |/i|. Our results extend the findings of Ref. [15]. A related expansion valid in the 
limit of an approximate U{1)pq symmetry can be found in Ref. [7]. 
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The CP-odd Higgs bosons of the NMSSM are made up of Im(S)/\/2 and the non- 
Goldstone combination A° of Im{H^)V2 and Im{H^)/V2. In the {Im{S)/V2, A'^J basis 
the mass matrix is 

2 / ( ^ ] (4 + ah^ - 36 



f) (4 + «)7^ 



Mi = {^] ^''^ ■ ■ ,.„^w, I- (A.i: 




~ 


A 




XXv 


A = 




7 = 







K 




K /i 



In writing this expression, we have defined the variables 

XAx XA^ 

a = , d = . (A. 2) 

K fJ, K /i 

We will assume that 7, a, and S are all much less than unity, and treat A and S2/3 as being 
on the order of unityjj Under these assumptions, the ■M'a^^ element is much larger than the 
M.\^^ and ■M'aj^2 elements. Thus, the rotation angle to diagonalize the mass matrix will be 
small and the mass eigenstates will be close to Im{S)/\/2 and A^. Labelling these mass 
eigenstates by {a^, A"}, we find the mixing matrix to be 

2o ) - O ( M^yv^ ) . (A.3) 

with the approximate rotation matrix given by 

011 = l-4|(l-3a)7' = 022 (A.4) 

2A^ 

O,, = ^7f2-3a + 3a^-3^5V— 7' = -02i. (A.5) 

2A V A ; 2A3 

In deriving these expressions we have treated a, 7 = 0(e) and 5 = 0{e^) with e <^ 1, and 
we have kept terms only up to 0{e^). 

The corresponding mass eigenvalues to this level of approximation are 

2 

= -35 + 9a^7^ (A.6) 

2A 

^(l + «)+4^(l-2«)72. (A.7) 

S2f3 2A 

To discuss the NMSSM CP-even Higgs mass matrices and their eigenvalues, it is con- 
venient to work in the basis {/i°, h^,H'^}, where /i° is the combination of Re{H^)/^/2 and 
Re{H^) / V2 that has the same tree-level couplings to the electroweak gauge bosons as the 




This expansion also works quite well for small k/X and S2/3 up to higher-order terms in a. 
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SM Higgs, H^ is the orthogonal combination, and h'^ = Re{S) / \/2. The transformation to 
this basis is simply 



Re{Hl)/V2 \ 
ReW)/V2 
Re{S)/V2 J 



cos/3 sin/3 

— sin P cos P 

1 




(A. 



The CP-even symmetric mass matrix in the basis {h^, /i°, iJ^} reads 

^ 1\4b + 9^ cL) 2A7 - 7(2 + a)s2/3 7^(1 - 9^)s2l3C2f3 



M' 



H 



\ 



A + 5 + ^a^' 



V 



-7(2 + a)c2(3 



.(A.9) 



Here, we have defined 



f 



{9' + 9") 
2A2 



(A.IO) 



The {/i°, /i°, iJ°} basis is useful because all the mixing elements in this matrix are suppressed 
by at least a factor of 7, while the lower two diagonal elements are of order unity. Since this 
mixing is small, we will designate the mass eigenvalues by {/i°, hs-, H^} in analogy with the 
corresponding MSSM states. 



The transformation to the mass eigenbasis is 



h' \ 


. 


^ K 


hs 


= u 1 


K 


H' 


' 


K H?, 



with the unitary matrix U given by 



(A.ii: 
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f/ll 

f/l2 



r 

32 



2A - (2 + a)s2p 



u. 



21 



^((-4 + 5)[2A-(2 + a)s2^] 

3 

+Y^((A-S2/3)(3p-6As2/3 + sy 

-24^[2As2;3 + ^'(A - 2s2;3)(A + S2/3)] 
A 



-^((-4 + 5)([2A-(2 + a)s2^]) 
3 



2r2 

^'-2/3 



(A - 2s2pf 



2s2piy - 3As2;3 + s') + ~g'\{\ - ?>S2fi){\ - 2s 



2/3 J 



(A.12) 
(A.13) 



(A.14) 



f/' 
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1 - o7M (A - S2,3) [A - (1 + a)s2/3] 



■[(-l + a)A + 2(l + a)s2/3] 



4c2 p2 

^''2/3''2/3 



(A - 2S2/3)' 



(A.15) 



[/ 



13 



7 S2f3C2p 

4 A 



-2 + a)X + 2{a + g^ - a^ )s2/3 



U: 
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7^ / C2;3si^ 



2 VA(A-2S2;3) 



[A + ^^A-2^s 



2/3j 



-c^ - -'^'''' , [(2 + ~g')\' - 4f\s2, - 4(1 - ~g')s' 



U23 = 7 



A(A - 2S2;3)^ 

C2/3S2/3 

2(A-2s2^)2 

1 



2/3J 



[(2 - a)\ - (4 + 2a - 5)s2/3] + a' 



S2I3C2I3>^ 

2(A - 2s2^)3 



(A.16) 
(A.17) 



(A + 2S2;3) 



(A.18) 



KA - 2s2py 



-l \,, ^ ^, C2pS2p A^ - 2(4 + f)\h2p + (19 + 10^^)A2s^^ - 4(5 + 3^^)As^^ + 12s^^ 



U- 
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-7 1 ../^"!^" A {^-a)\-{A + 2a-8)s2,]+a' ''''''\ (X + 2s2,) 



+r 



2(A-2s2/3)2 
2A(A - 2s2,3)3 



U- 



33 



1+7' 



„2 2 

■^2/3^2/3 



2(A - 2s2^)3 



2(A - 2S2;3)3 

3c2^A - (A - 2s2;3)[(l + f)X' - 4^'As2/3 + 2fs%] 



(-l + a)A + 2(l + a)s2/3 
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(A.19) 



(A.20) 



The CP-even Higgs mass eigenvalues are 



mlo 



^Is 



m]jo 






g^c\p - A^ + (2 + q;)s2/3A - as 



-2/3 

4 + 5 

2 i \2 , '^■^2/3 



2/3 



(A.21) 
(A.22) 



/i 



+Y y + 



2A(l + a) 

S2/3 



2A 



(2 + a)As2/3 + (1 + a)s2fl - -^ 



2ci^g2/3 

(A - 252/3)- 



:[A-2(l + a)s2/3] 



(A.23) 



2C2^S2;3 



+f I (-1 + ~g')sl^ + -_l^p— [A - 2(1 + a)s2;3] | • 

[A — 2S2/3J 



B RG Equations with Exotics 

We collect here the modifications to the RG equations that arise when multiple sets of vector- 
like exotics are added to the NMSSM. Our notation conventions follow NMHDECAY. The 
additional exotics we consider consist of D © D'^, and L © L'^, with D = (3, 1, —1/3) and 
L = (1, 2, —1/2). Taken together, these exotics have the quantum numbers of a 5 © 5 of 
SU{5). We include the trilinear exotic superpotential couplings 



WD^D,SD,D': + ^^^SL,t.. 



We also add the soft terms 



^soft 



X^nfi C mjj^\Di\' + mj^c\D^i\' + ml\Lj\' + mic\L^ 



c|2 



HdAd^ SD, Dt + ^l,Al^ SL, ^. 



(B.l) 



(B.2) 
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The RG equations for the superpotential couphngs are 
dlnX 



(IGtt^ 



dt 



'IQn' 



dluK, 
dt 



:i67r^) ^^ 
' ^ dt 



-2 [{Yl + Yl)f + 2C./] 

3iV,A2 + 6«:2 + 3iV, ^ ^d, + 3iV^ ^ ^1,), 



AT^A^ + 2«:2 _ 2 [(F^^ + F^.)/ + 2 C,gl 



'167r' 



dt 



N^\^ + 2«:2 - 2 (F^ + y2 )^/2 ^ 3 ^^(^ 



2 A„/2 



+24+^c5:el,, + iV.E4- 



Here A^c = 3 is the number of colors and iV^ = 2 is the number of "weak" colors. 
For the trilinear A terms, we have 
dAx 



Wtt' 



dt 



(167r 



'IQn^ 



2\ dAf^ 

dAp^ 
dt 



2{N^ + 2)X^A^ + 4k^A^ + 2N,{X^tAt + A^Ab) + 2A^A, 
+4 [{Yl + Y^^)g'^M, + 2C2g^M^] 

+2N, Y, ilAD. + 2iV^ Y. 4^^- 
« i 



2N^\^A^ + 4«:2a, + 4 (F^ + F^ )(?'2^i + 2 Cg^^g'M; 



+4^1,, Ab, + 2K J2 ^IM + 2iV- E 4^^. 



'IGtt 



dt 



2N^yA^ + 4/s:^A, + 4 {Y' + F^^c)(7'^Mi + 2 C^g'M^ 



-^H.Al, + 2iV, 5^ HAd. + 2iV^ 5^ 4 Al, . 



(B.3) 
(B.4) 

(B.5) 

(B.6) 



(B.7) 
(B.8) 

(B.9) 

:b.io) 
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Finally, the modifications to the running of the scalar soft masses is 



dt 



+2N, Y, i% i^l. + ^% + ^1 + ^ A ) 

i 

+2iV^ ^ i%{m% + ml. + m| + ^i^), (B.li: 



3 



(JTX) 

(levr)^— ^ = 2e^^(m^^+m^. + m| + A^J (B.12) 



dt 



^{Ykg"Mi + C3<73'M3^) + 2F^X, 



(167r)2^ = 24K^, + mi. + m| + Ai; (B.13) 

The quantity C in these expressions is the hypercharge Z)-term which vanishes in many simple 
models of gauge mediation [18j . 
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